Local Superfluidity of Parahydrogen Clusters 
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We study by Quantum Monte Carlo simulations the local superfluid response of small (up to 
27 molecules) parahydrogen clusters, down to temperatures as low as 0.05 K. We show that at 
low temperature superfluidity is not confined at the surface of the clusters, as recently claimed by 
Khairallah et al. [Phys. Rev. Lett. 98, 183401 (2007)]. Rather, even clusters with a pronounced 
shell structure are essentially uniformly superfluid. Superfluidity occurs as a result of long exchange 
cycles involving all molecules. 

PACS numbers: 67.90.+Z, 61.25.Em 



Small clusters of parahydrogen (P-H2) molecules are 
believed to remain liquidlike at significantly lower tem- 
perature than that at which bulk P-H2 crystallizes, ren- 
dering the observation of superfluidity (SF) possible 
P, [3,13, 01 ■ Indeed, some experimental evidence of SF 
has been obtained for clusters of approximately 15 
molecules, at a temperature r=0.15 K Q . Those findings 
prompted an ongoing theoretical research effort, aimed 
at characterizing structural and superfluid properties of 
g-H? clusters, bo th p ristine and doped with impurities 
E a i, S, 0, S, S, [ifl El, [13 ■ The theoretical method of 
choice has been computer simulations based on Quantum 
Monte Carlo (QMC) techniques, at zero or finite temper- 
ature, which yield essentially exact numerical estimates 
for Bose systems. 

Although numerical discrepancies exist among the var- 
ious works, attributable in some cases to the use of dif- 
ferent interaction potentials, a finite global superfluid re- 
sponse has been generally predicted for pure P-H2 clus- 
ters of size at least up to N=27 molecules, at T < 1 
K; some clusters (of size 22 < < 30) turn superfluid 
as they undergo quantum melting, going from solid- to 
liquidlike when T is lowered, purely as a result of quan- 
tum exchanges; the solid bulk phase is found to emerge 
non-monotonically for increasing cluster size [3, [3] ■ 

One interesting open issue is how the superfluid re- 
sponse is distributed across the droplet. Because a cluster 
is a flnite non-homogeneous system, one may meaning- 
fully pose the question of whether the decoupling from an 
externally applied rotation characterizing a superfluid is 
greater in some parts thereof (e.g., the surface). The rel- 
evance of this aspect goes obviously beyond the physics 
of P-H2 clusters, as it impacts the general theoretical un- 
derstanding of SF, and its connection with Bose conden- 
sation. 

In a recent paper. Khairallah et al. utilized the radial 
density proflle of particles involved in permutation cycles 
to study the local superfluid density of P-II2 droplets. 
Based on their results, they concluded that the superfluid 
response is largely confined at the surface of these small 
clusters, mainly arising from exchange cycles involving 
loosely bound surface molecules [3|- 



Aside from the difficulty of providing a meaningful def- 
inition of "surface" of clusters of such small size, the 
physical picture proposed in Ref. [3l is at odds with 
recent simulation work by us, on pure and isotopically 
doped P-H2 clusters, which has yielded indirect evidence 
that their superfluid response correlates with permuta- 
tional exchange cycles compri sing molecules in both the 
inner and the outer shell [J[3, [ll[. Indeed, even clusters 
that display a coexistence of liquid- and solidlike behav- 
ior at T 1 K [e.g., (^-112)23] are essentially entirely 
superfluid in the liquidlike "phase" [3, [3| , observation 
which raises doubts about the contention that SF may 
be confined at the surface of the cluster. 

In order to provide an independent theoretical check 
of the predictions of Ref. [3] , we have undertaken in this 
work a systematic study of the local superfluid properties 
of P-H2 clusters at low temperature, using QMC simula- 
tions with a recently introduced, rigorous local super- 
fluid density estimator, consistent with the basic two- 
fluid model [l3| . Specifically we have computed the local 
superfluid response of pristine P-H2 clusters of up to 27 
molecules in the temperature range 0.05 K < T < 3 K, 
and examined its connection to global properties previ- 
ously discussed in Refs. 0, [3]- 

Our results show that, contrary to what stated in Ref. 
[1], SF in these clusters is not localized at the surface; 
rather, in the low temperature limit clusters are uni- 
formly superfluid; this statement applies not only to the 
more liquidlike ones (i.e., those with fewer than 20 
molecules), but also to those displaying a more pro- 
nounced shell structure, which softens at low tempera- 
ture due to quantum melting, in turn underlain by ex- 
change cycles involving molecules in both the first and 
second shells. Based on these findings we contend that 
the physical picture proposed in Ref. [3| of SF as arising 
from "loosely bound surface molecules" does not apply. 
Indeed, we have observed no clusters with a rigid core 
and a superfluid outer shell. Instead, SF crucially de- 
pends on the onset of long exchange cycles involving all 
molecules, not just those on the surface. 

We adopted the usual microscopic model for our sys- 
tem of interest, namely a collection of A^ P-H2 molecules 
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regarded as Bose particles of spin S=0, interacting via 
a central pair potential. We used the most commonly 
adopted model interaction for P-H2 clusters, namely the 
Silvera-Goldman potential , in all the calculations for 
which results are shown here. We studied this system 
by QMC simulations at finite temperature based on the 
continuous-space worm algorithm 3, 1^ j which gives us 
access to temperatures as low as 0.05 K; technical details 
of our calculations are the same of Ref. 0]. Profiles of 
radial superfluid density were computed using the local 
estimator proposed in Rcf. [l^; specifically, the local 
superfluid response ps{r) is expressed as follows: 

ps{r) Am'^{AA{r)) 



p{r) 



(1) 



where m is the mass of a P-H2 molecule, j3 = 1/T, p{r) 
is the local density and (■ • ■) stands for the thermal av- 
erage. A is the total area swept by the many-particle 
paths projected onto a plane perpendicular to one of the 
three equivalent rotation axes, whereas A(r) and Icif) 
are, respectively, the contributions to the total area A 
and to the classical moment of inertia Ic from a spheri- 
cal shell of radius r centered at the center of mass of the 
cluster. The above estimator is obtained from the def- 
inition of normal fraction as proportional to the linear 
response of a system to an externally applied rotation. 
The response of the system as a whole is decomposed in 
separate contributions, which provide local information. 

Figure [1] shows profiles of total [p{r)\ and superfluid 
[p5(r)] radial densities, computed with respect to the 
center of mass of the system, for a cluster comprising 
A^=18 P-H2 molecules. At T = 1 K the cluster is entirely 
superfluid, within the statistical uncertainty of our calcu- 
lation [3, Q . As shown in Fig. [H the superfluid fraction 
is 100% everywhere in the cluster, with no sign of weak- 
ening near the center. At T = 2 K, the total superfluid 
fraction of the cluster decreases to ~ 55%; notably, how- 
ever, the superfluid density remains finite throughout the 
cluster, even in its inner region. Upon increasing T to 3 
K, the superfluid fraction drops to about 14% and ps(r) 
is correspondingly uniformly depressed throughout the 
whole system, almost completely in correspondence of 
the total density minimum. The latter information con- 
firms the qualitative relation between SF and exchange 
cycles involving molecules in different spatial regions of 
the cluster. Clearly, the local superfluid response of this 
cluster is never conflned to the surface. It is important to 
note that the discussed variations of ps{r) occur against 
a total density profile which is essentially unaffected by 
temperature, for this liquidlike cluster (as shown in Fig. 

HD. 

A behavior analogous to that illustrated above for the 
(p-H2)i8 cluster is displayed essentially by all clusters 
with fewer than 22 molecules, which are liquidlikc in the 
temperature range considered here. As discussed in Rcfs. 
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FIG. 1: (color online) Profiles of total [p(r)] and superfluid 
[pg(r)] density computed with respect to the center of mass 
of the cluster (p-H2)i8 at T = 1 K (lower panel), T = 2 K 
(middle panel) and T — 3 K (upper panel). 



0.04 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 

r(A) 

FIG. 2: (color online) Profiles of total and superfiuid density 
at T = 0.0625 K (solid line and circles), T = 0.25 K (dashed 
line and triangles) and T —0.5 K (dotted line and diamonds) 
for the cluster (p-H2)26. 

0, 3 , clusters with a number of molecules between 22 
and 30 feature a different physical behavior, displaying 
(to different degrees) coexistence of liquid- and solidlike 
properties at T < 1 K. We discuss here the specific case 

(P-H2)26- 

Figure [2] shows the same profiles as Fig. [T] but for a 
cluster of 7V=26 molecules. This cluster displays insu- 
lating properties at T=0.5 K; its superfluid fraction is 
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barely over 10%, rising to approximately 32% at T = 
0.25 K and reaching 100% at T 0.0625 K. Correspond- 
ingly, and quite unlike the N=18 case, the structure of 
the cluster undergoes a significant change, as shown by 
the evolution of the radial density profile. In particular, 
the peak corresponding to the first shell (at r ~ 2 A) be- 
comes depressed as T decreases; hardly any remnant of 
the shell structure is left at T = 0.0625 K, as a liquidlike 
structure emerges, i.e., the cluster undergoes quantum 
melting 0, [^. Indeed, it is at low T (in contrast with 
what suggested in Ref. [1]), when this cluster quantum 
melts, that its basic physics becomes qualitatively remi- 
niscent of that of a helium cluster. 
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FIG. 3: (color online) Frequency of occurrence of permutation 
cycles involving M molecules in a cluster of 26 P-H2 molecules 
at T=0.5 K (lower values) and T=0.0625 K (higher values). 
Long permutation cycles set in at low T. 



It is interesting to examine the corresponding evolu- 
tion of the local superfluid density, also shown in Fig. 
[2j At the highest temperature shown, ps(r) is small at 
all distances from the center of mass, and negligible be- 
tween the two shells (diamonds); here too it has to be 
stressed how, though the superfluid response is fraction- 
ally greater in the immediate surface (i.e., r > 7 A), 
it is nonetheless not confined thereto. Moreover, as the 
temperature is lowered and the cluster turns liquidlike 
and superfluid, the relative increase of ps(r) is far more 
significant in its inner part. Especially telling is the dis- 
appearance of the dip in the region between the two shells 
(at r 4 A), suggesting the occurrence of permutations 
involving particles in both the inner and outer shell. This 
is consistent with the observed temperature dependence 
of the frequency of occurrence of permutation cycles in- 
volving different numbers of molecules, shown in Fig. [31 

We discussed in previous work the non-monotonic 
behavior of the total superfluid fraction of clusters of size 
A'^ > 22; we further illustrate this phenomenon here, by 
focusing on a very narrow range of cluster sizes, within 
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FIG. 4: (color online) Profiles of total and superfiuid density 
for the cluster (p-H2)25 (solid line and circles), (p-H2)26 (dot- 
ted line and diamonds), (P-H2) 27 (dashed line and triangles) 
at T = 0.5 K. Open symbols show estimates of the local su- 
perfluid density obtained using the same measure utilized in 
Ref. 0. 



which the value of the superfluid fraction changes dra- 
matically on adding just one molecule at a time. The 
results for the local superfluid, and total density defini- 
tively demonstrate that such behavior rcfiects structural 
changes throughout the whole cluster, not just its surface 
as contended in Ref. [1]. 

Total and superfluid densities for clusters of 25, 26 and 
27 molecules at a temperature of 0.5 K are shown in Fig. 
m The superfluid fraction is close to 1 for A^ = 25, drop- 
ping to roughly 0.12 for A^=26 and rising again up to 
approximately 0.85 for N=27. For the N=25 cluster, the 
superfiuid response is fractionally close to 100% through- 
out the entire cluster. While obviously fewer particles 
reside on average in the core than on the surface region 
of such a small cluster, the results of Fig. 3] again clearly 
show how SF is not confined to the surface, and the core 
of the cluster is not rigid but superfiuid. The addition of a 
single molecule causes the entire cluster to turn abruptly 
insulating and solidlike, as shown by the increase in the 
height of the total density first shell peak. The local su- 
perfluid response is concurrently depressed throughout 
the whole system, especially in the region between the 
two shells, to indicate the paucity of exchanges taking 
place between molecules in different shells. On adding 
yet another molecule, i.e., on increasing the cluster size 
to N=27, the cluster returns to a liquidlike structure, and 
consistently the superfiuid response is enhanced, at the 
surface, but also in the inner part of the system. These 
findings lend substantial theoretical support to the no- 
tion that the liquidlike, superfiuid phase of these small 
clusters occurs as a result of the onset of long exchanges 
involving all molecules in the system. This is in funda- 
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mental disaccord with the picture put forth in Ref. [j], 
according to which SF would correlate to long exchanges 
involving loosely bound surface molecules. 

One might be tempted to attribute the discrepancy 
of our results with those of Ref. Q to the difference 
between the estimator of the local superfluid density uti- 
lized here, and the qualitative measure employed in Ref. 

which lacks any microscopic justification 13, T^]. For 
comparison, we show in Fig. [4] a superfluid density pro- 
file obtained using the same measure of Ref. Q (open 
circles); no substantial difference is seen with respect to 
our result (full circles). Furthermore, our physical con- 
clusions can also be inferred from the evolution of the 
density profiles as the temperature is lowered (Fig. [2]), 
clearly showing the disappearance of the shell structure 
and the concurrent approach to unity of the total su- 
perfluid fraction (for this quantity the same estimator is 
used in our work and Ref. 4), as well as by the temper- 
ature dependence of the frequency of occurrence of long 
permutation cycles (Fig. [3]). All of this unambiguously 
shows that superfluidity correlates with long exchanges, 
comprising molecules both in the inner and outer shells. 

More generally, it should be pointed out that the 
continuous-space Worm Algorithm, utilized in this work 
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a much greater efficiency in sampling per- 
mutations than the method used in Ref. 0. It is quite 
conceivable that the physical behavior observed in this 
work, underlain by long permutation cycles, may have 
remained hidden to the authors of Ref. 0, due to the 
inefficiency of their method in sampling such cycles, es- 
pecially at the lowest temperatures. 

In order to explore a possible dependence of the phys- 
ical results on the potential utilized, wc have also per- 
formed in this work calculations with the Buck poten- 
tial [3] (another popular choice), as well as with the 
Lennard- Jones potential utilized in Ref. [J. Our main 
physical conclusions are unaffected by the intcrmolccular 
potential utilized, which essentially only influences the 
energetics and the temperature scale. 

Summarizing we have studied the size and tempera- 
ture dependence of the local superfluid density of P-H2 
clusters. We found no evidence of superfluid response 
confined at the surface, regardless of temperature and 
system size explored. SF concomitantly grows as T de- 
creases in the inner and outer part of the system in liq- 
uidlike clusters [i.e., (p-Il2)i8] as well as in those which 
feature quantum melting at low T. In the latter case, 
structural variations such as the progressive collapse of 
the first peak of the radial density profile are due to long 



(i.e., inter-shell) exchanges, involving all molecules in the 
cluster. Long exchanges cause the system to turn liq- 
uidlike at sufficiently low T, with the ensuing onset of 
superfluidity. Significant supcrfiuid response is only seen 
when the the core and the surface are both superfluid and 
liquidlike, i.e., it is not the effect of quantum exchanges 
involving exclusively, or mostly, surface molecules. 
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